Elastic heterogeneity of soft random solids 
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Spatial heterogeneity in the elastic properties of soft random solids is investigated via a two- 
pronged approach. First, a nonlocal phenomenological model for the elastic free energy is exam- 
ined. This features a quenched random kernel, which induces randomness in the residual stress 
and Lame coefficients. Second, a semi-microscopic model network is explored using replica statisti- 
cal mechanics. The Goldstone fluctuations of the semi-microscopic model are shown to reproduce 
the phenomenological model, and via this correspondence the statistical properties of the residual 
stress and Lame coefficients are inferred. Correlations involving the residual stress are found to be 
long-ranged and governed by a universal parameter that also gives the mean shear modulus. 

PACS numbers: 61.43.-j,62.20.Dc,82.70.Gg 



Introduction — As a consequence of randomness incorpo- 
rated at synthesis, random solids (e.g. polymer networks, 
glasses, a-Si) are heterogeneous. For example, the mean 
positions of the constituent particles exhibit no apparent 
long-range order, and every particle inhabits a unique 
spatial environment. Particularly for soft random solids, 
such as rubber, in which the particle positions undergo 
large thermal fluctuations, heterogeneity also manifests 
itself via the RMS particle-displacements, which are ran- 
dom and continuously distributed 0, • 

The elasticity of rubber, and especially its softness 
with respect to shear deformations, have been studied for 
many years via the classical theory, developed by Kuhn, 
Flory, Wall, Treloar and others Q and based on a mi- 
croscopic picture of Gaussian polymer chains. While the 
classical theory has proven highly successful, it is a homo- 
geneous theory, preserving no information about the ran- 
dom structure of these essentially heterogeneous solids. 
Thus, it is incapable of describing consequences of the 
heterogeneity, such as random spatial variations in the 
local elastic parameters and residual stress, or non-affine 
deformations that occur in reponse to applied stresses Q . 

The purpose of the present Letter is to present a theo- 
retical development that goes beyond the classical theory 
of rubber elasticity by accounting for the heterogeneity. 
What emerges is an elasticity theory featuring spatially 
fluctuating Lame coefficients and residual stresses, to- 
gether with a statistical characterization in terms of their 
mean values and spatial self- and cross-correlations. In 
particular, we find that not only is the stress-stress cor- 
relation long ranged — behavior that can be argued for 
on general grounds — but so are the cross-correlation be- 
tween the residual stress and, e.g., the shear modulus. By 
contrast, we find the self-correlation of the shear modu- 
lus to be short-ranged. Furthermore, we find that the 
long-ranged correlations and the average shear modulus 
are governed by a common universal parameter that is 
independent of microscopic details. 

To obtain our statistical characterization of soft ran- 
dom solids we take the following route. First, we ex- 



amine a nonlocal phenomenological model of a random 
elastic medium, which we subsequently derive from a 
semi-microscopic model. We then determine the state 
to which it relaxes when randomness is present, and re- 
expand the elastic free energy around this new equilib- 
rium reference state This relaxed state is, however, 
still randomly stressed Q; nevertheless, the stress in the 
relaxed state — the so-called residual stress — satisfies the 
mechanical equilibrium condition djO~jk{x) = 0. In its 
local limit, the proposed phenomenological model re- 
produces a version of Lagrangian elasticity theory that 
features random Lame coefficients and residual stresses. 
Second, we consider the statistical mechanics of a cer- 
tain semi-microscopic model of a random-solid-forming 
system — the randomly linked particle model (RLPM) y| • 
Via the replica method, applied to the RLPM to deal 
with its structural randomness, followed by an analysis 
of Goldstone fluctuations of the random solid state of the 
RLPM, we arrive at precisely the aforementioned nonlo- 
cal phenomenological model, except that the latter has 
been appropriately disorder-averaged using replicas. By 
comparing these two disorder-averaged models, we char- 
acterize the elastic heterogeneity of soft random solids. 
Phenomenological model — We begin by examining the 
following nonlocal model for the elastic free energy T 
of soft random solids, which describes "mass points" 
(i.e. coarse-grained volume-elements) that interact with 
one another through random harmonic attractions: 

r = \ J dz\ dz 2 G(zi, z 2 )^\R(z 1 ) - i?(z 2 )| 2 - \zi - j 
+ ±A [ dz (\dR(z)\ - if , (1) 



i.e., a functional of the deformation of the system R(z), 
which specifies the D-dimensional position vector to 
which the mass point at position z is displaced. 
denotes the determinant of the deformation gradient ten- 
sor dRi/dzj and, correspondingly, the parameter Ao, 
which we take to be large, heavily penalizes density vari- 
ations. It results from a competition between (i) re- 
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pulsions (either direct or mediated via solvent, e.g., 
excluded- volume) and (ii) intcrmolccular attractions and 
external pressure. The nonlocal kernel G(zi,z 2 ) de- 
scribes link-induced harmonic attractions between mass 
points, originating in the entropy of the molecular chains 
of the heterogeneous network, and are modeled as a 
"zero-rest- length" spring of random spring coefficient. 
G(zx, Z2) is taken to be a quenched random function of 
the two positions, z\ and Z2, symmetric under z\ <-> z 2 . 
This model, which we shall shortly see emerging from the 
RLPM, is in the spirit of the classical theory of rubber 
elasticity 

The free energy V is not stable at the state R(z) = z for 
two reasons: first, the attraction G causes a small, spa- 
tially uniform contraction [the fractional volume change 
being O(l/Ao)]; second, the randomness of G addition- 
ally destabilizes this contracted state, causing the adop- 
tion of a randomly deformed equilibrium state. We de- 
note this relaxation as z — > z = C,z + v(z), in which 
£ describes the uniform contraction, and v(z) describes 
the random local deformation. This process can be un- 
derstood in the setting of the preparation of a sample of 
rubber via instantaneous crosslinking: crosslinking not 
only drives the liquid-to-random-solid transition but also 
generates a uniform inward pressure, as well as intro- 
ducing random stresses. As a result, immediately af- 
ter crosslinking the state is not stable, but relaxes into 
a new equilibrium state determined by the particular 
realization of randomness created by the crosslinking. 
Free energy stationarity applied to Eq. (JJJ shows that 
at large Ao the contraction £ « 1 — (p/DXa), where 
P = h I dz2 ( Zl ~ Z2 ) 2 G o{zi - z 2 ), and G (zi - z 2 ) = 
[G(zi,z 2 )] is the disorder average (denoted by [■■•]) of 
G, and we have assumed that this disorder averaged G 
is translationally and rotationally invariant. As we shall 
see below, p is actually the mean shear modulus. 

The random relaxation v(z) is, in general, difficult to 
determine, owing to the large nonlincarity from the Ao 
term. However, by assuming that G\ = G — Go (the 
random part of G) is small, and seeking v(z) to first 
order, free energy minimization ST/Sv(z) = enforces 
that the Fourier transform v(p) = J dz exp(-ipz) v(z) is 
given by 



v(p) 



P h (p) • f(p) 



r T (p)-f(p) 

2(Go(0)-G (p)) X P 2 + 2(Go(0)-G (p)) ' 



(2) 



where f(z) = — 2£ J dz\ G%(z, Zx) (z — z%), and V^ k {p) = 

Pj Pk/P 2 and 'PjkiP) = S jk ~ PjkiP) are ' respectively, 
transverse and longitudinal projector, and j, k, etc. are 
Cartesian indices. 

Next, we transform to a new reference state, which 
we take to be the relaxed state, via the coordinate 
transformation z — » z{z) = Qz + v(z) [with inverse 
z{z)\. Correspondingly, we obtain the relaxed-state ker- 
nel G(zi,z 2 ) = G(z(z{), z(z 2 )). We then re-expand the 
elastic free energy as a functional of the deformation 



R(z) = z + u(z) around the new reference state, keep- 
ing terms to quadratic order in the nonlinear Lagrangian 
strain tensor £jk(z) = \{{dRi/ dzj){dRi/ dzk) — Sjk) and 
to sub- leading order in 1/Ao- By considering the local 
limit of the resulting free energy via a gradient expan- 
sion, and dropping an additive constant, we arrive at 

r=y > rfz{Tr ( j(5)-e(z^ M (5)Tr£(5) 2 +^(Tre(5)) 2 }. (3) 

This free energy features three random elements: the 
heterogeneous Lame coefficients p(z) and A(z) and the 
residual stress ajk(z); respectively, they are given by 

p(p) = pS(p)-iC 1 p- 2 (p-f(P)), (4a) 
X(p) = X $(p) + 2(iC 1 p- 2 (p ■ f(p)) -pi®), ( 4b ) 



cr-jk 



(P) 



G\(p — q,q) + iSjk 



P ■ f (P) 

CP 



n2 



+*c 1 r 2 {p 3 vl (p) + p k vjM)fx (P), (4c) 

where 3{p) = (27r)- D 5(p). It can be shown that a obeys the 
equilibrium condition dj(Tjk(z) — 0. (Strictly speaking, 
the residual stress a, defined here, is the leading-order 
term of the true stress, in the sense of a gradient expan- 
sion, owing to the non-locality of our model. As a result, 
the equilibrium condition only holds to the corresponding 
order in the gradient expansion.) 

To determine the statistics of G, we first treat the 
statistical mechanics of elastic fluctuations governed by 
Eq. using replicas. [We explain why this procedure is 
appropriate, below, before Eq. (| 1 4|> . ] Thus, we arrive at 
the disorder-averaged replicated partition function 



i zn ] = /ni =1 pircx p(- r »/ r )' 



(5) 



where we have set k-Q — 1 and a labels replicas; [Z n ] 
is related to the disorder averaged Gibbs free energy via 
[T] = -Tlim„^ ^M^"] W- The resulting effective 
pure free energy is 



P AO 

1 n — 9 



\dR a \ 



\\ [G(z 1 ,z 2 )]^ R (z 1 ,z 2 ) 



s , , [G(zi, z 2 )G(z 3 , z 4 )] c * R (z 1 , z 2 )^r(z 3 , z 4 )H ,(6) 

'{*«} 

where the subscripts {z,i} indicate integration variables, 
'c' indicates a cumulant (or connected correlator), and 



9 a (z u z 2 ) ee {\R a (zi) - R a (z 2 )\ 2 - \zx - z 2 \ 2 }. (7) 

a=l 

The dots represent terms arising from higher-order cu- 
mulants, which characterize the non-Gaussian nature of 
the distribution of G. 

Semi-microscopic model — Our next goal is to determine 
the cumulants via the statistical-mechanical analysis of 
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a semi-microscopic model. Thus, we consider the ran- 
domly linked particle model (RLPM) Q, which consists 
of N particles having coordinates {cj}^ =1 interacting via 
an excluded-volume term, all in a fluctuating volume 
V, the mean value of which is controlled by a pressure 
p. In addition, permanent soft links, modeled as "zero- 
rest-length" springs, are introduced at random, with a 
separation-dependent probability, between nearby parti- 
cles in a liquid-state configuration, so that the number 
of links is quasi-Poisson-distributed and the correlations 
among the links are consistent with the correlations of 
the liquid state. This is a version of the Deam-Edwards 
distribution [Toj , and is at least appropriate for instanta- 
neous schemes for synthesizing soft random solids. 

For a given realization of the randomness, the Hamil- 
tonian of the RLPM is given by 



H 



RLPM 



* — ^i.-j — l * — *e= 1 



M K 



2a 2 



(8) 



The first term describes excluded-volume repulsion [y 
is taken to be large and, thus, density variations are sup- 
pressed); the second term describes attractions associ- 
ated with the M randomly-chosen pairs {ci c , Cj e }^Li of 
particles that have been softly linked to one another. 

The particles of the RLPM can be identified with small 
molecules or coarse-grained polymers, and the soft links, 
e.g., with Gaussian molecular chains. In making this 
coarse-graining one is assuming that microscopic details 
(e.g., the precise positions of the crosslinks on a polymer, 
the internal conformational degrees of freedom of the 
polymers, and the effects of entanglement) do not play 
significant roles for the long wave-length physics. In part, 
these assumptions are justified by studying more detailed 
models, in which the conformational degrees of freedom 
of the polymers are retained Q- Thus, the RLPM can 
be regarded as a caricature of vulcanized rubber or as a 
model of chemical gels or other soft random solids. It 
is a model very much in the spirit of lattice percolation, 
except that it naturally allows for particle motion as well 
as connectivity, and is therefore suitable for the study of 
continuum elasticity and other issues associated with the 
(thermal or deformational) motion of particles. 

To analyze the RLPM we use replica field theory, which 
allows us to average over the quenched randomness, and, 
accordingly, we introduce the following order parameter, 
suitable for describing the soft random solid state: 



a))] 



Here, x = {x°,x x ,...,x n ) denotes the (1 + n)-fold- 
replicated D-vector x, and (■ • • ) denotes a thermal av- 
erage for a given realization of randomness (i.e. a given 
configuration of the links) 0- This order parameter is 
capable of detecting and characterizing the random solid 
state. The additional (i.e. th ) replica arises from the 
preparation ensemble associated with the Deam-Edwards 
distribution of the quenched randomness; this replica has 



a fixed volume Vq. We fix the pressure in the measure- 
ment ensemble to equal the mean pressure in the prepa- 
ration ensemble, i.e., p = v ^ r i + ^rf-- This leads to the 

ZV Q Vq 

following effective pure ^-dependent free energy: 



n+n 



NTrj 2 
2V n A 



]T A(p)n(p)n(-p) 



pGHRS 

AT In / ,/,-exp ( 



nv 2 N 2 
2V 



npV 



— Y 

yn 

pGHRS 



A(p)n(p)e^ £ ). (9) 



Here, £l(p) and A(p) are, respectively, the Fourier trans- 
forms of Q(x) and A(x) (= exp(— x 2 /2a 2 )), and A(p) = 
n" =0 A(p Q ) and A = A(p)| p=0 . Then, the disorder- 
averaged Gibbs free energy is given by 



[Q] 



d r 

-T lim — — In / Dft exp (—Qi+ n 

n^O dn J 



IT). 



(10) 



The parameter r/ 2 controls the density of links, with the 
average number of links to any single particle being r/ 2 . 
The restriction p e HRS serves to exclude macroscopic 
density fluctuations, i.e., to account for the effects of the 
strong excluded-volume interactions @. At the level of 
mean-field theory, the transition from the liquid to the 
random solid state occurs at rj 2 = 1, as can readily be 
seen via an expansion of Qi+n in powers of fl. For larger 
rj 2 , the equilibrium value of the order parameter is 

n(*)=0|||drp(r)e-«^-«l"-E^I« B -M"} > (ll) 

where z = (z, Qz, . . . , Qz) represents the random mean 
positions about which replicated particles undergo ther- 
mal fluctuations, ( — (V/Vq) 1 ^ is the contraction ratio, 
and Q and p(r) characterize the state via, respectively, 
the fraction of localized particles and the distribution 
of (inverse square) localization lengths. Again within 
mean-field theory, Q, p and £ are determined by self- 
consistency conditions that follow from making Gi+n< sta- 
tionary; in particular, Q obeys 1 — Q = exp(— n 2 Q) UQj- 
With the approach to the RLPM in place, we now de- 
scribe how the equilibrium value of the order parameter is 
modified when the system undergoes an elastic deforma- 
tion. In this replica theory, such deformations amount to 
a replica version of Goldstone excitations. In view of the 
pattern of spatial symmetry breaking, they are specified 
by n (not n + 1) independent long wavelength deforma- 
tion fields {R 1 (z) , R 2 (z) , . . . , R n (z)\, each depending on 
a single replica of space z [ill Il2j , with the constraint 
det (dR a /d((z)) = 1. This constraint corresponds to the 
low-energy excitations of the model: pure shear deforma- 
tions from the contracted measurement state. Hence, we 
arrive at the appropriately deformed order parameter: 



-)exp(-r|x-i?(z)| 2 /2), (12) 



where R(z) = (z, R 1 (z), . . . , R n (z)) are the random mean 
positions of the particles in the preparation (a = 0) and 
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(deformed) measurement (a = 1, . . . , n) replicas, and x 2 
denotes ^2™=o\ xa \ 2 - The um<; density and excluded- 
volume parameters will be taken to obey: v 2 N/TVq 3> 
rj 2 » 1. We take rj 2 > 1 (and hence Q m 1) be- 
cause we shall be concerned with the well-linked regime, 
in which only a very small fraction of particles are not 
part of the infinite cluster. As we shall see, the con- 
dition v 2 N/TVq ^> 1 indicates that the medium is near- 
incompressible, and v 2 N/TVq 3> ij 2 holds because we are 
concerned with soft random solids, for which the bulk 
modulus greatly exceeds the shear modulus. 

To determine the free energy of the Goldstone deforma- 
tion we insert the deformed order parameter into Gi+n, 
thus obtaining (up to an irrelevant additive constant) 

v 2 N 2 (V V , f 

Ql+n[R] = ~eyrr [ TT ~ 1 + \ \ K^Zi, Z 2 ) * r{z U Z 2 ) 
~WT K 2 (zi,Z 2 ,Z 3 ,Z i )^ R (z ll Z 2 )^ R (z 3 ,Z i )^ .(13) 

J{*} 

The kernels K%(zi, z 2 ) and K 2 [z\, z 2 , 23, Z4) are given by 
rather lengthy formula; in terms of rj 2 , Q and p, but 
are, in essence, bell-shaped functions of the separations 
of the variables that fall off on the scale of the typical 
localization length; we shall report on them elsewhere. 
Results — Observe that Eqs. 10 and (|13|) are equivalent, 
in the sense that they are both capable of describing the 
free energy (disorder-averaged via replicas) of spatially 
varying deformations and homogeneous volume varia- 
tions in soft random solids, about an unrelaxed reference 
state R{z) = z. This equivalence makes the RLPM a sen- 
sible candidate for addressing the question posed above 
about the statistics of G\. By comparing the these two 
equations we arrive at the correspondence 

[G{z!,z 2 )] = Kiizuzz), (14a) 
[G(z 1 ,z 2 )G(zs,Zi)] c = K 2 (z 1 ,z 2 ,z 3 ,z 4 ), (14b) 
A = u 2 N 2 /V 2 , (14c) 

Higher-order cumulants of G may be determined via the 
expansion of G[R] to higher order in "3/^, as we shall 
undertake elsewhere. 
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TABLE I: Long- wavelength variances and covariances of 
the elastic properties of soft random solids in the relaxed 
state. The entry in row R and column C, when multiplied 
by (N/Vo)T 2 d(p + p'), yields the connected disorder correla- 
tor [R(p) C(p% = [R(p) C(p')] - [R(p)] [C(p% 

Now that we have ascertained information about the 
statistics of G, we use it to address the statistics of the 
heterogeneous Lame coefficients p and A and the residual 
stress a of the relaxed state, which we do via Eqs. J2J) and 



(JIJ. This lead to the following results for the disorder- 
averaged elastic parameters: 

{[a jk {z% [p{z% [\{z)]) = (0,p,v 2 N 2 /V 2 ), (15) 

where the average shear modulus p is given by |l3j 

p = 6TN/V , 6 = -±7 1 2 Q 2 + e- r ' 2 Q + r ] 2 Q~l. (16) 

(rj 2 — l) 3 near the transition, and ~ ?/ 2 for rj 2 3> 1. 
The variances and covariances among the Fourier trans- 
forms of p,, A and a at long wavelengths are given in 
Table ||| In these formulae, 7 is a scalar determined by 
rj 2 and p; the tensor A is a p-dependent structure that 
vanishes on contraction with p. 

By transforming these variances and covariances back 
to real space, one can obtain their leading large-distance 
behavior. It is evident that the entities not involving the 
stress a{r) (i.e. [p(0)p(r)} c , [A(0)A(r)] c and [/i(0)A(r)] c ) 
are short-ranged in real space: more precisely, they are 
proportional to 5 s (r) (i.e. to the Dirac delta function 
smoothed on the scale of the short-distance cutoff, which 
should be taken to be the typical scale for localization 
length, in order to validate the Goldstone-fluctuation 
framework). By contrast, those entities involving the 
stress have more interesting behavior: in three dimen- 
sions and at large-scales we find that 

[<M0)<7«(r)] c = ^T 2 ^B ljkl /r 3 , (17a) 
[a ii (0) M (r)] c =-^8T 2 ^ (V^(r) - |^(r))/r 3 , (17b) 
M0)A(r)] c = ±8T 2 % (?>L (r) _ i-p?.(r))/r 3 , (17c) 

where the tensor is a complicated structure com- 

prising terms built from projection operators of the vec- 
tor r, together with various index combinations. It is per- 
haps unexpected that the correlation between the resid- 
ual stress and shear modulus is a long-ranged quantity. 

By starting with a microscopic model, we have con- 
structed a random nonlocal elasticity theory, together 
with a statistical characterization of the parameters that 
feature in this theory. So far, we have focused on this 
statistical characterization in the local limit of the the- 
ory, Eq. and have found — inter alia — that a uni- 
versal parameter 9 controls all long-range correlations. 
This parameter also governs the large-distance statistics 
of non-affine deformations that occur in response to ap- 
plied stresses |6J . The statistical content and implications 
of the theory can also be explored beyond the local limit. 
For example, the complete statistics of the nonlocal ker- 
nel G are amenable to the present formalism, and will 
allow the extension of the statistics of non-affine defor- 
mations down to the mesh scale. 
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